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ABSTRACT 

We compare the Li abundances of a sample of stars with planets discovered with the 
Dopplcr method to a sample of stars without detected planets. We prepared the sam- 
ples by combining the Li abundances reported in several recent studies in a consistent 
way. Our results confirm recent claims that the Li abundances of stars with planets 
are smaller than those of stars without planets near the solar temperature. We also 
find that the vsini and i? H x anomalies correlate with the Li abundance anomalies. 
These results suggest that planet formation processes have altered the rotation and Li 
abundances of stars that host Doppler detected planets. We encourage others to test 
these findings with additional observations of Li in stars with temperatures between 
5600 and 6200 K. 



1 INTRODUCTION 



Stellar Li abundances are at once very informative and dif- 
ficult to interpret. This follows from the relative delicacy of 
Li nuclei in the shallow surface layers of stars, where they 
are destroyed via (p, a) reactions when they are mixed to re- 
gions with warm protons. Li abundances in dwarf star photo- 
spheres are observed to correlate with effective temperature 
(T e ff), age, rotation, binarity and metallicity. However, even 
within a single coeval population (such as the open cluster 
M67) variations in Li abundance are observed among stars 
that otherwise appear identical ( Randich et al.|2006 |. 

Some have suggested that the presence of a protoplan- 
etary disk is the missing parameter that accounts for the 
observed spread in Li abundances among similar stars. The 
the surface Li abundance of a star could be altered via ac- 
cretion of protoplanetary disk material (Gonzalez 19981 or 



a significant depletion of Li among SWPs relative to a com- 
parison star sample, but only in the T e ff range 5600 to 5850 
K. |Takeda fc Kawanomoto] ( |2005| ) largely confirmed their 
finding s for the T eft range 5800 to 5900 K. |Chen fe Zhao] 
( |2006[ ), restricting their attention t o T eff = 5600 to 590 
K, also confirmed the conclusions of |Israelian et aL] (2004) 



However, Luck & Heiter (20061, employing a larger com- 



parison star sample, did not find a significant difference be- 
tween SWPs and a comparison sample. They attribute the 



Li abundance difference found by Israelian et al. (20041 to 
a systematic difference in the temperature scales in their 
study and the study of |Chen et a l. (2001), the results of 
which they had used to supplement their small comparison 
sample. Most recently, Takeda et al. (2007a I measured Li 



via a change in its rotation ( Chen & Zhao 2006 Takeda et 
al.|2007a |. Thus, while Li has the potential to serve as a use- 



ful probe of stellar and planetary processes, their effects on 
Li abundances are difficult to disentangle given our present 
level of understanding. 

Nevertheless, several studies have attempted to isolate 
the effects of planets on Li abundance. Gonza lez fc Laws] 
12000) first suggested that stars with planets (SWPs), when 
corrected for simple linear trends with T e g, metallicitjrland 



in 118 nearby solar analogs. While they included only a few 
SWPs in their study, they again concluded that SWPs tend 
to have smaller Li abundances. 



The purpose of the present study is to resolve the con- 
flicting findings concerning Li abundances in SWPs. It con- 
tinues our series of studies on the chemical abundances of 
nearby SWPs (for a summary of previous papers, see |Gon-| 
zalez & Laws (2007l). In Gonzalez & Laws (20071 we corn- 



age, display smaller Li abundances than field stars. Ryan 
( 2000 1 examined the Li abundance trends more carefully and 



concluded that any possible differences were not significant. 
|Gonzalez et al.| ( |2001[ ), employing a larger sample, agreed 
with his conclusion. 



Israelian et al. ( 2004 \ revisited this topic and reported 



1 In the following we use either [Fe/H] or [M/H] as a proxy for 
stellar metallicity. 



bined the chemical abundance data from multiple studies 
in a consistent way and compared the chemical abundances 
of SWPs and comparison stars for several elements; we did 
not include Li in the comparison, because it requires a fun- 
damentally different analysis. We employ similar methods 
in the present study. In Section 2 we present new samples 
of SWPs and comparison stars formed by combining the Li 
abundance results from multiple studies. In Section 3 we use 
these new samples to determine whether SWPs have differ- 
ent Li abundances than stars without detected planets. We 
also examine vsini and the i? H K activity index in SWPs. We 
discuss our findings within the context of planet formation 
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models in Section 4. We present our conclusions in Section 
5. 



2 COMBINING SAMPLES 



In Gonzalez & Laws (2007) we presented the results of our 



chemical abundance analyses of 18 elements (including Li) 
in 31 SWPs. We also calculated corrections needed to place 
our data and the data from other similar studies on the 
same abundance scale. This procedure allowed us to con- 
struct larger and more accurate SWP and comparison star 
samples, which we employed to compare C, O, Na, Mg, Al, 
Si, Ca, Sc, Ti and Ni abundances. We concluded that there is 
evidence for significant differences between the two samples 
for some elements. 

We adopted a similar procedure for Li. We begin with 
the Li abundance values listed in Table 1 of Gonzalez fc Lawsl 



(20071 for 31 SWPs. To these data we added Li abundance 



values (and upper limits) of other SWPs from our previous 
papers. We excluded SWPs with T cfi < 5000 K and log g 
< 4.0. The derived stellar properties tend to be less accu- 
rate for cooler stars, and few stars with T e g < 5000 K have 
detectable Li. The limit on surface gravity restricts the sam- 
ple to main sequence stars; a star that has evolved off the 
main sequence onto the subgiant or giant branch undergoes 
dramatic changes in its internal structure that results in 
changes in the surface Li abundance. Having applied these 
limits, the total number of SWPs with Li abundances (and 
upper limits) from our studies comes to 52. 

Next, we compiled Li abundances (and upper limits) for 
SWPs with T cff > 5000 K and log g > 4.0 from the follow- 
ing studies (number of SWPs retained from each study is 
also indicated) : Israelian et al. ( 2004 \ 
( 2006[ ) - 48, |Takeda fc Kawanomoto| ( |2005 
|et al.| ( [2007a| ) - 5. The same studies include Li abundances 
for stars without planets; the totals (after applying the same 
limits) are:|Israelian et al.|(|2004| - 34,|Luck fc Heiter| (|2006[) 
- 112,lTakeda fe Kawanomotof ( [2005 I - 98, |Takeda et ah] 
( 2007a| - 113P| We selected these studies for analysis be- 
cause they have similar temperature scales, and their data 
are of comparable quality. 

One complication that we did not have to consider for 



83, Luck & Heiter 



27, Takeda 



the elements examined in Gonzalez & Laws (20071 is the 



presence of upper limits in the abundance data. The mini- 
mum detectable Li abundance varies among the studies we 
employed, because their spectra have differing S/N ratios. 
The Li abundance values from our series of studies typically 
have the lowest upper limits, while those from |Luck fc Heiter| 
( 2006 1 have the highest (but the upper limits from the other 



studies are not much lower). We show the stars from |Luck fc| 
Heiter ( 2006 ) that have Li abundance upper limits in Figure 



1. We also include in the figure a straight line that matches 
the upper envelope of the upper limits for T e fl < 6200 K; 
it is described by the equation log Li = —3.41 + 0.00079167 
T c ff. We will employ this line as a cutoff for selecting Li 

2 Note, we removed HD 33636 from the list of SWPs, given the 
finding by |Bean et af] | |2007| l that its companion has a mass in 
the stellar regime. The SWP totals listed for each study include 
HD 33636, but we moved it to the stars without planets category 
prior to conducting the analysis in the next section. 
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Figure 1. Li abundance upper limits versus T e g for SWPs (dots) 
and stars without planets (open circles) from |Luck fe Heiter] 
pOOtit . The upper limit cutoff curve is shown as a dashed line. 



abundance values to include in our final sample preparation 
below. 



Next, we selected the Luck & Heiter ( 2006 1 study as the 



reference and adjusted the data from the other studies to be 
consistent with it. First, we adjusted the T e ff values in each 
study by applying a simple linear equation calibrated us- 



ing SWPs in common between it and|Luck fc Heiter ( 2006 1 



We then adjusted the log g values using a linear equation 
including a term with the adjusted T e rf values. The adjust- 
ments to [Fe/H] and Li are simple constant offsets; typical 
[Fe/H] and Li abundance offsets are about 0.02 and 0.05 
dex, respectively. 

The cutoff line shown in Figure 1 is only valid for 
T c ff < 6200 K, since two of the stars without planets from 



Luck & Heiter ( 2006 I above this temperature have upper 



limits much higher than predicted from the extrapolated 
cutoff line. One possible reason for this discrepancy is the 
fact that faster rotation is more common among earlier spec- 
tral type dwarfs, resulting in broader, shallower absorption 
lines. Therefore, in our second cut of the data, we removed 
stars with T c ff > 6250 K; we also excluded stars with 
T c ff < 5550 K, given the paucity of SWPs with detected 
Li in this range. Our upper T e ff limit also excludes stars 



that might display the Li dip phenomenon (see Balachan- 



dran ( 1995 l), which would further complicate our analysis 



below. Finally, we excluded all the stars with log Li values 
falling below the cutoff line (whether upper limits or detec- 
tions). 

For those stars present in multiple studies, we calculated 
simple averages of their parameters. We show the resulting 
Li abundances plotted against T c ff in Figure 2. The final 
samples include 37 SWPs and 147 stars lacking planets. 

The corrected individual T c ff values differ from the |Luckl 
|fc Heiter] ( |2006[ ) values with a typical dispersion of about ± 
70 K. The corresponding dispersion for log g is about ± 
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Ax = 30 ! log T^° mp - log T P ff | + |[Fe/H] comp - [Fe/H] p | 
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Figure 2. Final adopted Li abundances versus T e fj. Symbols 
have the same meanings as in Figure 1. 



0.13 dex, while for [Fe/H] and log Li they are ± 0.03 and ± 
0.08 dex, respectively. These numbers are consistent with the 
quoted typical measurement uncertainties in the individual 
studies. Of course, the standard error of the mean (s.e.m.) 
values for the average values of these quantities for those 
stars included in multiple studies are smaller. 



3 COMPARING SAMPLES 
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Two stars with identical values of T e fr, log g, [Fe/H] and 
M v will have a Ai value of zero. Therefore, two stars with 
identical ages, masses and compositions will also have a Ai 
value of zero. Our motivation in choosing these particular 
coefficients is to give comparable contributions to Ai from 
each term for two stars that differ by about 2a in each pa- 
rameter. We increased the contribution of the log T e fi term, 
however, to account for its greater importance in determin- 
ing the Li abundance. And, we reduced the contribution of 
log g, given the relatively lesser importance of this quantity. 

We calculated Ai for each SWP relative to each com- 
parison star, resulting in 5439 values. For each SWP, we then 
selected the comparison star with the smallest value of Ai. 
We plot the differences in Li abundances between each SWP 
and the closest matching comparison star in Figure 3a. 

One weakness of this approach is that only a small sub- 
set of the data is used to calculate the differences in the Li 
abundances. These results are sensitive to outliers. We can 
improve upon our analysis as follows. In our second approach 
we calculated weighted average values of Li abundance dif- 
ferences for each SWP using all the comparison star data. 
The weights are given by (Ai) -2 . We show the results of 
this analysis in Figure 3b. 

The low Li abundances evident for the SWPs near 5800 
K in Figure 3a are also evident in Figure 3b. The average Li 
abundance difference for SWPs with 5800 <T cff < 5950 K 
is —0.38 ± 0.10 (s.e.m.) dex. Figure 3b also shows a general 
excess of Li among SWPs with T G ff > 5950 K, which have a 
mean Li abundance excess of 0.12 ± 0.04 (s.e.m.) dex. 

To test the sensitivity of our results to the form of the 
A index, we repeated the above analysis with the following 
modified version: 



Objective comparisons between the Li abundances of SWPs 
and comparison stars have proven to be difficult. Some stud- 
ies have done nothing more than simple visual comparisons 
on a plot like Figure 2, while others have attempted to cor- 
rect for differences in T e ff, [Fe/H] and chromospheric ac- 
tivity (e.g., Gonzalez fc Laws] ( |2000[ >) . We propose a new 



approach to comparing these samples. 

Observations of stars in open clusters and in the field 
indicate that Li abundances are most sensitive to T e ff and 
age (see Cutispoto et al. (20031; Takeda et al. (2007a I and 



references therein). In addition, the rate of Li destruction in 
a star's envelope should be sensitive to [Fe/H], since a star's 
convection zone is deeper for larger [Fe/H] ( |Montalban &c 
Rebolo|2002 1. For each star in our samples we have T e ff, log 
g, [Fe/H] and M v , all corrected for small systematic differ- 
ences among the various studies. 

Isochrone-based stellar ages can be derived from the 
T e ff, [Fe/H] and M v values, and this has been done in sev- 
eral studies of SWPs (e.g., |Takeda et aL|p007b)). Uniform 



age estimates are not available for all the stars in our SWP 
and comparison star samples. However, our method of com- 
parison does not require that we calculate ages. Instead, we 
calculate an index that is a measure of the proximity of two 
stars in T e g, log g, [Fe/H], M v -space. It is based on the 
following equation: 



A 2 = 40 | log T^° mp - log T P ff | + |[Fe/H] comp - [Fe/H] p | 
+0.5 | logg comp - logg p | + 1.5 |M V ™ P - M P | 

We show the resulting Li abundance differences calcu- 
lated using the A2 index in Figure 4. There are only a few 
differences from Figure 3, but the pattern of low Li abun- 
dances among the cooler SWPs and the high Li abundances 
among the hotter SWPs evident in Figure 3b remain un- 
changed. The averages and dispersions are nearly identical 
to those calculated with the Ai index. 

Rotation is another important parameter that corre- 
lates with Li abundance (Cutispoto et al. 20031. Uniform 



rotational velocity (vsini) measurements are available from 
a single source for most of the stars we plotted in Figures 
3 and 4. Valenti & Fischer (2005) measured vsini for the 
1040 stars in their SPOCS survey; they quote a typical un- 
certainty of 0.5 km s _1 |^]This level of precision is adequate 
to resolve the range of vsini values measured for SWPs (< 1 
to ~ 15 km s -1 ). Of the 37 SWPs plotted in Figures 3 and 



3 |Valenti fc Fischer| ( |2005} assigned a value of km s — 1 to 70 
stars in their SPOCS database; 31 one of them pass our T e ff 
cuts. We assigned a value of 0.3 km s _1 to these stars prior to 
calculating the A indices. 
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Figure 3. Li abundance differences (SWPs - comparison stars) 
versus T e g. Panel a shows the differences between SWPs and 
the most similar comparison stars. Panel b shows the weighted 
average differences. See text for additional details. 
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Figure 4. Same as Figure 3 but using the A2 index. 



4, 33 have vsini values listed in Valenti & Fischer (20051; in 
addition, 70 of the 147 stars in our comparison sample have 
vsini values listed in their study. We define a new A-index 
that includes the vsini term: 



A 3 = 30 j log T™ mp - log T P ff | + |[Fe/H] comp - [Fe/H] p | 
+0.5 I logg comp - logg p | + \M c v ° mp - M P | 
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Figure 5. Same as Figure 3 but using the A3 index. 
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We show the Li abundance differences calculated us- 
ing the A3 index in Figure 5; the pattern apparent in 
the previous figures is still present. The average Li abun- 
dance difference for SWPs with 5800 < T cff < 5950 K is 
-0.31 ± 0.10 (s.e.m.) dex; SWPs with T cff > 5950 K have 
an average Li abundance difference of 0.15 ± 0.05 (s.e.m.) 
dex. 

We repeated the calculations using an index with the 
vsini term three times the weight of that in the A3 index. 
The resulting Li abundance differences are very similar to 
those in Figure 5. 

We summarize the properties of the SWP and compar- 
ison stars plotted in Figures 3 to 5 in Table 1. 



4 DISCUSSION 

We confirm recent claims that the Li abundances of SWPs 
with T e ff near 5800 K tend to be lower than those of stars 
without detected planets. We also find, for the first time, 
evidence that SWPs hotter than about 5900 K have excess 
Li. The magnitudes of the average the Li abundance differ- 
ences for the cool and hot SWPs are larger than the typical 
measurement errors. We cannot reach any firm conclusions 
regarding SWPs cooler than 5800 K, given the paucity of 
such stars with detectable Li in our sample. 

While our method of analysis does not rely on stellar 
evolutionary models, it does permit us to minimize the ef- 
fects of differences in age and mass in our analysis of the 
Li abundances of the SWP and comparison star samples. 
In addition, the trends apparent in Figures 3 and 4 per- 
sisted even after correcting for differences in vsini. Inclusion 
of vsini as a parameter reduced the magnitude of the aver- 
age Li abundance deficit among the cool SWPs only by 0.07 
dex. 
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Table 1. Properties of SWPs and comparison stars plotted in Figures 3 to 5. 
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Figure 6. Plot of vsini differences versus T e g for the same SWPs 
plotted in Figure 5. 
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Figure 7. Plot of vsini differences versus T e jj for SWPs and 
comparison stars from |VaIenti &: Fisc her ( 2005 ) . 



Although the effect is small, vsini apparently does have 
an effect on the Li abundances of SWPs. To further explore 
the relationship between vsini and the presence of planets, 
we calculated vsini differences using the Ai index. We show 
the results for the same SWPs plotted in Figure 5 in Figure 
6. Figure 6b shows that the vsini values are smaller than 
those of the comparison stars among the cooler stars; they 
are nearly the same for SWPs near 6000 K and much larger 
for the hottest ones. 

We can improve our comparison of vsini values by em- 
ploying the SWP and comparison star samples from |Valenti| 
& Fischer ( 2005 I . Applying the same T e gf and log g limits we 



used to prepare our samples above, we produced a samples 
of 82 SWPs and 596 comparison starsQWe show the result- 
ing vsini differences calculated in the same way as in Figure 
6 in Figure 7. These larger samples confirm our results from 
Figure 6. From the data in Figure 7b, we find that SWPs 



4 Since |Valenti fc Fischer| | |2005| was published, several stars orig- 
inally in their comparison stars list have been found to host plan- 
ets. We have moved such stars to the SWP sample. 



with T e ff between 5800 and 5950 K have an average vsini 
value of —0.77 ± 0.19 (s.e.m.) km s _1 relative to the com- 
parison stars; the corresponding average difference for stars 
hotter than 5950 K is +0.45 ± 0.39 (s.e.m.) km s" 1 . 

The pattern of Li abundance differences between SWPs 
and comparison stars shows a remarkable correspondence 
with the pattern of differences in vsini. Both Li and vsini are 
smaller for cool SWPs and larger for hot SWPs. These re- 



sults confirm the preliminary finding of Takeda et al. ( 2007a I 



that solar analogs with slower rotation have lower Li abun- 
dances. They reached this conclusion after having conducted 
detailed spectroscopic analyses of 118 solar analogs. They 
measured line widths, which allowed them to determine 
vsini+macroturbulence for each star. Unlike Valcnti & Fis- 
cher (2005), however, they did not determine separate vsini 



and macroturbulence velocities. What's more, they only in- 
cluded 5 SWPs in their sample, limiting their conclusions 
about the possible relationship between SWPs, rotation and 
Li abundance. 

Stellar chromospheric activity is known to correlate 
with rotation. If the pattern in vsini values we found among 
the SWPs is a real effect, then there should be similar pat- 
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Figure 8. Plot of log -R HK differences ve rsus T e ff for SWPs and 
comparison stars from Valcnti & Fischer ( 2005 ) and Wright ct al. 
( |2004| ). 
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Figure 9. Same as Figure 8 but excluding stars with logi? HK > 
-4.43. 



terns in measures of chromospheric activity. We can test 
this prediction using the -R H k chromospheric activity in- 
dex, which has been measured for many Sun-l ike stars. The 
best s ource of -Rhk values for our purposes is |Wright et al.| 
(20041. They tabulate -Rhk values for the stars included in 



the Valenti & Fischer ( 2005 1 study. In our first cut of the 



data, we cross referenced the stars in Wright et al. (20041 
with those in | Valenti fc F ischer (2005) to produce new sam- 
ples of SWPs and comparison stars. After applying our T c ff 
and log g cuts, we produced samples of 52 SWPs and 411 
comparison stars. We show the resulting R H k differences in 
Figure 8, calculated using the Ai index. As expected from 
the vsini results, SWPs tend to have smaller -Rhk values. 
From the results shown in Figure 8b, we calculate that the 
average difference in logi? HK for SWPs with T e ff between 
5800 and 5950 K is -0.10 ± 0.02 (s.e.m.) dex. 

We should note that some scatter in the vsini and -Rhk 
difference plots is expected. First, the equatorial velocity of 
a star is viewed from random angles. Second, like the Sun 
most FGK dwarfs display variable -R H k- 

Before we can attribute these results to physical causes, 
we need to consider the possible effects of biases. One im- 
portant bias is the variation in the S /N ratio of the spectra, 
which we already discussed in Section 2. We dealt with it by 
excluding stars with Li abundances below the level where 
they could be uniformly determined. 

Other biases are caused by differences in vsini and chro- 
mospheric activity, which affect the detectability of planets 
with the Doppler method (all the SWPs included in our com- 
parisons above were discovered with the Doppler method). 
A large vsini value results in broad lines, and a large value of 
-Rhk is accompanied by greater photospheric "jitter," which 
is a source of noise in Doppler measurements (Wright 2005). 
The selection criteria of planet search groups includes a re- 



quirement that the spectra be sufficiently sharp-lined to de- 
tect the small Doppler shifts induced by an orbiting planet. 
Spectroscopic analyses of the type used to derive the ba- 
sic stellar parameters (including Li abundance) also require 
sharp-lined spectra. Therefore, both the SWP and compar- 
ison star samples exclude stars with large vsini and R H k 
values. 

In order to test for possible SWP selection bias resulting 
from stellar activity, we produced a second set of SWP and 
comparison stars samples. This time we excluded stars with 
log R HK > —4.43 dex; this is the largest value of logR H K 



measured for an SWP (HP 22049;|Jenkins et aL] f2006) quote 
a value of —4.43 dex, while Wright et al. ( 2004 1 quote a value 
of —4.51 dex). This second set contains the same number 
of SWPs (HD 22049 is too cool to have been included in 
our SWP samples) and 393 comparison stars. We plot these 
data in Figure 9. The average difference in logR H K f° r the 
SWP data in Figure 9b with T off between 5800 and 5950 
K is —0.07 ± 0.01 (s.e.m.) dex. Therefore, the activity bias 
can only account for a small part of the difference in R H k 
between SWPs and comparison stars. 

It is interesting to note that the star with the largest 
positive anomaly in Figures 8b and 9b, HD 179949, is the 
best candidate for star-planet magnetic interactions. [Shkol- 
nik et al. ( 2003 1 reported significant synchronous enhance- 



ment of Ca II H and K emission with the planet in orbit 
around the star. 

We also consider the vsini bias to have negligible effect 



on our results. First, the stars in the Valenti & Fischer ( 2005 1 
SPOCS database have median and maximum vsini values of 
2.4 and 54.7 km s _ , respectively, and only 44 stars have 
vsini values greater than 11 km s _ . Only about 5 per cent 
of the comparison stars have vsini values greater than that 
of the SWP in the SPOCS sample with the largest vsini 
value (HD 120136 with vsini = 15 km s _1 ). In addition, all 
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the stars in the SPOCS database with vsini values above 
16.5 km s" 1 are beyond the upper cutoff in T e ff (6250 K) 
that we employed in our sample preparations. 

Second, given that it would be more difficult to de- 
tect planets around stars with the highest vsini values, we 
would expect the SWPs sample to have a smaller average 
vsini value than a comparison sample that included stars 
with higher vsini values. However, we actually find that 
SWPs have positive vsini anomalies for T fr > 6100 K, the 
very range having the highest vsini values in the SPOCS 
database. 

Third, vsini tends to decline with age. Therefore, if the 
SWP sample is significantly older than the comparison star 
sample, then the observed vsini values of the SWPs will 
be smaller than those in the comparison sample. We have 
greatly minimized this possible source of bias in our analysis 
above. By weighting the vsini and Li abundance differences 
according to the similarities of the stars, we are, in effect, 
comparing stars that are similar in age. 

|Valenti fc Fischer] ( |2005[ ) also excluded spectroscopic 
binaries and chromospherically very active stars. Again, we 
don't consider these biases to significantly affect our results 
for the reasons we gave above. 

Are the Li abundance differences caused by differences 
in rotation or are they both caused by something else? The 
results shown in Figure 5b imply that vsini cannot fully 
explain the differences in Li abundance. The more likely ex- 
planation for the correlation between Li and rotation among 
SWPs is that both are the result of process related to planet 
formation. 

This is not the first time Li abundance has been linked 
with rotation in stars. The low Li abundances observed over 
a narrow temperature range centered near spectral type F5 
on the main sequence (called the Li dip) are accompanied by 
low vsini values; a similar phenomenon is observed among 
giants and subgiants at lower T c g in the Hertzsprung gap. 



Bohm-Vitense (20041 suggests that the simultaneous de- 



crease of both Li and rotation velocity in these two classes 
of stars is caused by evolutionary changes leading to deep 
mixing, which destroys Li at the higher temperatures at the 
base of the envelope and redistributes angular momentum. 
In addition, Li is known to be preserved in short period 
synchronously rotating binaries, confirming stellar models 



that include rotationally induced mixing (Ryan & Deliyan- 
|nis|1995[ ). 

Clearly, processes that affect mixing at the base of con- 
vection zone must be included in any theory offered to ex- 
plain the correlation between Li abundance and rotation. 



Israelian et al. (2004) proposed two hypotheses to account 



for the low Li abundances they measured for SWPs with 
T c ff near 5800 K. First, they suggested that since the proto- 
planetary disk around a young star contains a large fraction 
of the system's angular momentum, it can cause rotational 
breaking of the star. The breaking, in turn, results in deeper 
mixing in the star and thus more efficient destruction of Li. 
Second, in a system with a migrating giant planet, tidal 
forces from the planet could create a shear instability at the 
interface between the convective and radiative zones, leading 
to more mixing and hence destruction of Li. 

Our results for SWPs near 5800 K confirm the correla- 
tion between Li abundance and rotation based on the first 
hypothesis. Testing the second hypothesis would require ex- 



amining the orbits of the planets as well. To be effective, 
such an analysis should employ a larger sample than ours. 

Any theory must also account for the higher Li abun- 
dances and faster rotation we found for the SWPs near 6200 
K. In this temperature regime the mass of the convection 
zone is much smaller. Such a star might have weaker cou- 
pling to the protoplanetary disk. In addition, a less massive 
convection zone would be more susceptible to alteration of 
the surface composition by accretion of refractory material 
(|Gonzalez|p98t. 



Laws & Gonzalez (2003) and Ashwell et al. (20051 ar- 



gued that the anomalous chemical abundance pattern (in- 
cluding extremely high Li and Be abundances) of the early 
F dwarf J37 in NGC 6633 is best explained by accretion 
of chondritic material. This "self-enrichment" process could 
account for the Li excesses (though much smaller in magni- 
tude) we found for the late F SWPs. In this context it is no- 
table that Santos et al. ( 2004 ) found Be abundance excesses 
near 0.1 dex relative to comparison stars for SWPs between 
6000 and 6200 K. This is comparable to the Li excesses we 
determined for the same stars. Unfortunately, they only ob- 
served two comparison stars in this temperature range, so 
their results should be considered as preliminary. Finally, 
accretion should also increase the 6 Li abundance in a star's 
atmosphere, possibly to detectable levels ( Reddy et al.|200"2 
Israelian et al.|2003 1 



The self-enrichment hypothesis could also account for 
the higher vsini values of the hotter SWPs. A migrating 
planet orbits its host star in the same direction as it ro- 
tates. As material falls onto the star (whether small plane- 
tary building blocks or planets), it adds angular momentum 
to the star's envelope. 



5 CONCLUSIONS 

We prepared new SWP and comparison star samples by 
combining published spectroscopic data in a consistent way. 
Our analyses confirm those recent studies that have reported 
smaller Li abundances for SWPs near 5800 K. In addition, 
the data suggest that SWPs with T e ff near 6100 K have ex- 
cess Li abundances. The transition occurs near T c ff = 5950 
K. 

We also find that the patterns of vsini and -Rhk anoma- 
lies with T c ff correlate with the Li anomalies; the vsini and 
i? HK values are smaller among the cooler SWPs. 

The low Li abundance, vsini and R aK values among the 
cool SWPs is consistent with the suggestion that a proto- 
planetary disk causes rotational breaking of its host star, 
leading to additional mixing in its envelope. The additional 
mixing, in turn, accelerates the destruction of Li in the 
star's envelope. On the other hand, the preliminary find- 
ing of high Li abundance and vsini values of the hot SWPs 
can be counted as evidence for the self-enrichment hypoth- 
esis. Additional measurements of Be and 6 Li in SWPs and 
stars without planets near T e ff = 6100 K can help us test 
the self-enrichment hypothesis for SWPs in this temperature 
range. 

The number of SWPs with detected Li is still relatively 
small. The present situation can be improved by conducting 
detailed spectroscopic analyses on the more recently discov- 
ered SWPs. In addition, several published studies of SWPs 
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are based on spectra of modest S/N ratio (~ 150), prevent- 
ing detection of Li in many of the cooler stars. It would be 
straightforward to obtain spectra of higher S/N 300) for 
these stars. It would also be helpful to observe additional 
stars without detected planets with comparable S/N ratio. 

Our findings bring to five the number of stellar param- 
eters that correlate with the presence of Doppler detected 
planets: metallicity, mass, Li abundance, vsini and -R H k- 
There might also be additional secondary parameters, such 
as the 6 Li/ 7 Li isotope ratio and Al/Fe and Si/Fe abundances 
ratios, but these require confirmation. Taken together, these 
five parameters should make it possible to select stars with 
a high probability of hosting planets. 
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